Laser Induced Selective Alignment of Water Spin Isomers 
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We consider laser alignment of ortho and para spin isomers of water molecules by using strong 
and short off-resonance laser pulses. A single pulse is found to create a distinct transient alignment 
and antialignment of the isomeric species. We suggest selective alignment of one isomeric species 
(leaving the other species randomly aligned) by a pair of two laser pulses. 



PACS numbers: 33.80.-b, 37.10.Vz, 42.65.Re 
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I. INTRODUCTION 



According to quantum mechanics, a molecule that con- 
tains two identical atoms whose nuclei have nonzero spins 
can exist in the form of spin isomers. In particular, a wa- 
ter molecule exists in one of the two spin isomers, ortho 
or para with parallel or antiparallel proton spins, respec- 
tively. The ortho-para conversion is highly improbable in 
the gas phase [1] . Water isomers separation (or their ratio 
measurement) can be of great importance in medicine, bi- 
ology, chemistry, magnetic resonance imaging, etc. Mea- 
surement of the para to ortho ratio was used to estimate 
the formation temperature of molecules in space A 
method for mixture enrichment in one of the water iso- 
mers was reported using selective adsorption onto various 
surfaces Q, although some of these results were hardly 
reproducible. Recently Andreev et al [3] discussed an ex- 
pected difference in the response of the ortho and para 
forms of water to a non- uniform dc electric field near solid 
surfaces, which may potentially lead to isomers separa- 
tion. We are currently exploring an alternative route to 
ortho and para separation based on selective field-free 
alignment and antialignment of water isomers by short 
nonresonant laser pulses (for a review on field-free align- 
ment see Q). By alignment we mean the angular lo- 
calization of a molecular axis along the laser polariza- 
tion direction, while antialignment assumes localization 
of the same molecular axis in the plane perpendicular to 
the laser polarization. Although previous experimental 
studies on dissociative fragmentation by relatively long 
laser pulses have not found alignment of water molecules 
0, we show that a sizable alignment can be achieved in 
ultra-short regime. The selective alignment was recently 
experimentally demonstrated in a mixture of ortho and 
para modifications of linear molecules (^^A'^2 considered 
as an example) 0. Here, for the first time to our knowl- 
edge, we consider selective alignment for bent tri-atomic 
molecules like water. The paper is organized as follows. 
We briefly discuss the classification of the states of the or- 
tho and para water molecules in Sec. [ill The interaction 
of the water molecules with the laser field is considered 
in Sec. IIIII The results are presented in Sec. IIVI and 
briefly summarized in Sec. |Vl In the Appendix, we give 
a detailed description of the rotational wavefunctions of 
the ortho and para water isomers. 



II. CLASSIFICATION OF THE EIGENSTATES 
OF THE WATER MOLECULE 

The water molecule wavefunction is given by: 



\^^sp^n)\'^^ 



(1) 



where 



spill 

}, \^vib} ana \^eiei 
clear spin, vibrational and electronic wavefunctions, re- 
spectively. The water molecule eigenstates are classified 
using the C2v{M) Molecular Symmetry (MS) group 
i.e. using the irreducible representations Ai, A2, Bi, B2. 
We provide the group operations as well as its character 
table in the Appendix. According to the Pauli Exclu- 
sion Principle, the total wavefunction should change its 
sign after permutation of the hydrogen nuclei. Since we 
consider \^vib) and \^eie) in the ground state (which is 
symmetric with respect to the C2u(M) operations), only 
{^rot) and \^spin) are effected by the operations. The 
rotational part of the water molecule Hamiltonian is pre- 
sented by the rigid rotor model: 



and l^'eie) are rotational, nu- 
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where J is the angular momentum operator, (a, 6, c) are 
the molecule's principal axes, and la, h, Ic are the corre- 
sponding moments of inertia { la < h < Ic)- Rotational 
states of Al and A2 symmetries do not change sign as a 
result of permutation and belong to the para states (with 
anti-parallel spins), while the states of Bi and B2 sym- 
metries correspond to the ortho states. Further informa- 
tion about the water molecule eigenstates classification 
is given in the Appendix. 



III. INTERACTION OF THE WATER 
MOLECULES WITH LASER FIELD 

We consider the alignment process based on interaction 
with short off-resonance laser pulses that induce molecu- 
lar polarization, interact with it, and excite rotational 
wavepackets. The interaction of the water molecule 
(shown in Fig. [1]) with the laser pulse is given by Q: 



Hint = -^^£p{t)ap^p'ep'{t) , 



(3) 
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where the p index denotes the space-fixed Cartesian co- 
ordinates, Sp is the laser pulse envelope component, and 
a is the polariz ability tensor. 
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FIG. 2: The coordinate system. (X, Y, Z) are the space fixed 
laboratory axes, (a, b, c) are the body fixed principal axes, and 
[9, 4>, X) s^rs the Euler angles which define the body orientation 
in space with respect to the space fixed coordinates 



FIG. 1: The water molecule geometry. The axes are the prin- 
cipal axes of the molecule, a,b and c (a and b are in the 
molecular plane, and c is perpendicular to it). The prin- 
cipal moments of inertia are: la ~ 1.025 ■ lO'^^fcg ■ m^, 
h = 1.921 ■ 10"*^% ■ and h = 2.946 ■ 10"^^fcg ■ (note 
that Ic = la + lb according to ) 

For linearly polarized fields oriented in the laboratory 
Z direction, the interaction term is reduced to: 
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The Euler angles 9 and x, as well as the principal axes 
of the molecule (a, b, c) appear in Fig. [2] that presents 
the coordinates used in this work. In Eq. B, a"'' = 
c^aa — ctbb, etc, -0,*^^. are the rotational matrices and 

R={9,(b,x)- 

In the case of a symmetric top molecule (i.e. when 
ctxx = ctyy = a± and azz = ct\\), the Hamiltonian is re- 
duced to the well known form: Hind = ^ j^"^ (t) '^ct cos^ 9, 
where Aa — Q!|| — a±. In the case of water, the 

polarizability tensor components are: aaa = 1.528A , 
abb = 1.468A^ and a^c = 1.415A^ 

Given an initial rotational eigenstate of the asymmetric 
water molecule, we are interested in the field-free evolu- 
tion of the rotational subsystem after an excitation by 
a laser pulse that is short compared to the typical rota- 
tional periods of the water molecule. Using the impulsive 



approximation, we consider only the contribution from 
the interaction term to the Schrodinger equation during 
the pulse: 

{liD^oo + l2[Dl2+Dl_,]}^ = ih^ , (5) 

where 71 and 72 are found from Eq. ([¥]). 

Solving Eq. (O, we obtain the relation between the 
states before and after the pulse: 

ii>{t = 0+) = 

exp!^-^{l3,Dlo + f32[Dl,+Dl_^]}'^^{t = Q-) , (6) 

where /3i,2 = (l/h) J ji,2dt. 

In order to find ^'(t = 0+), we introduce an artificial 
parameter ^ that will be assigned the value ^=1 at the 
end of the calculations: 



exp|i{/3ii?2^+/32[i?^2 + ^o-2]}e}*(t = 0-) . (7) 

We express "i!^ in terms of the water molecule eigen- 
states: 



(8) 
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where J is the angular momentum quantum number, and 
m is the quantum number of the angular momentum pro- 
jected on the laboratory Z axis. For the chosen laser field 
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polarization, m is a conserved quantum number. This 
allows us to consider the excitation of the states with a 
different initial m-value separately. Index r (which is as- 
signed the values — J, ...,+J) numerates various 2 J + 1 
eigenstates arranged from the lowest energy (r = — J) to 
the highest energy (t = +J). 

Then we take a derivative of both sides of Eq. ^ with 
respect to ^, apply the bra (J',t',to| from the left, and 
obtain a set of differential equations: 



+P2{J'T'm\Dl2\jTm)+f]2{J'T'm\Dl_^\jTm)} , (9) 

where c = dc/d^. 

Since ^-^=0 = = 0") and ^-^^i = -^{t = 0+) (Eq. 
(Ill) ) , we solve numerically this set of equations from ^ = 
to ^ = 1, and find ^{t = 0+). 

The post-pulse field-free propagation of the created ro- 
tational wave packet is calculated as 



In order to calculate the alignment factor (cos^ d) at 
temperature T, one has to perform thermal averaging: 
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Here Qrot is the rotational partition function, ks is the 
Bolzmann's constant, {)t denotes quantum thermal aver- 
age over the molecular ensemble, and the indexes J, t, m 
denote the initial rotational state | J, t, m) of a molecule 
in the ensemble. The quantum average is calculated us- 
ing: 



{cos'' 0)jrUt) = (*(i)|cos2 0|*(i)> , 



(14) 



where 5'(<) is found by solving Eqs. (O, ^ and (fTO]) 
with the initial state ^'(t = 0^) = \J,T,m). The right 
hand side in Eq. can be easily calculated using the 
relation: 



(15) 
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where E'^'^ is the eigenvalue corresponding to the eigen- 
state I J, r, m) . We calculate the matrix elements in Eq. 
dl]) by expressing the water molecule eigenstates using 
the basis of the symmetric rotor: | J, k, m) ( J and m are 
defined as before, and k is the quantum number of the an- 
gular momentum projected on the rotor symmetry axis), 
and using the following relations: 



(Jkm\DlJj'k'm) ^ (-1)^=+™ X 

vwTwmT) I n (i I o . (11) 



Here s = 0, ±2, and ( "'^ "'^ "'^ | is the 3- j symbol, 

\mi 7712 ™3 / 

related to the Clebsch-Gordan coefficient by: 



Jl J2 J3 \ (-l)ii-J2-"3(2j + l)-l/2 X 
TOl TO2 TO3 J 

{jimij2m2\j3 - ms) . (12) 

The non-zero matrix elements from Eq. (fTT|) satisfy the 
following selection rules: |J— 2| < J' < J + 2 and k' = 
k,k±2. 

In order to evaluate the values of /3i,2, we consider the 
following pulse shape: e^{t) ~ eQexp|^^|. Consider- 
ing the maximal pulse intensity to be / = 3 • lO^^W/cm? 
and duration a = 20/s we evaluate Eq as Eq = 1.5 • 
10^°V/m and obtain /3i = -0.174 and (32 = -0.065. 



It is natural to consider the alignment of the para and 
ortho water molecules separately, since the laser pulse 
does not mix the para and ortho rotational states, i.e. if 
one considers an initial para rotational state and apply 
the pulse, only para rotational states are excited. The 
same holds for ortho rotational states. This is because 
the interaction term. Hind, in Eq- includes: 

DIq{R) cx (i?|J = 2,A: = 0,m = 0) 
^02(-R) + ^o-2(^) « (i?|2,2,0> + (i?|2,-2,0>,(16) 

which (according to Eq. in the Appendix) have Ai 
symmetry. In order to have a transition from a para 
to ortho state, the following matrix element should not 
necessarily vanish: {Para\Hind\Ortho) . Since the inter- 
action term has a symmetry Tind = ^1, the para state 
has Tpara — ^1 Or A2 symmetries, and the ortho has 
^ ortho — Bi or B2, direct product of the symmetries 
Tpara X ^ind X Tortho ucver coutains Ai and therefore nec- 
essarily {Para\Vind\Ortho) = (forbidden transition). 
Thus, the transitions between para and ortho states are 
forbidden and the sum in Eq. (|13p includes either initial 
para states or initial ortho states separately. 



IV. RESULTS 

Applying a linearly polarized 3-lO^^W/cm' , 20/s laser 
pulse, we expect the molecular highest polarizability axis 
(i.e. the a axis in Fig. [T|) to be aligned along the laser 
field polarization shortly after the pulse. Since the pulse 
is short compared to the typical rotational periods of the 
water molecule, it is considered as a delta-pulse. The cal- 
culated alignment factor (cos^ 9)t (Eq. is plotted 
in Fig. [3] as a function of time at temperature of 20K. 
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FIG. 3: The alignment factor after the laser pulse at tem- 
perature of 20K. Different dynamics of the ortho and para 
isomers is revealed during free propagation of the rotational 
wavepackets. 



FIG. 4: Alignment factor after the second pulse versus time 
at 20K. While the transient alignment of the para molecules 
is increased, the ortho alignment is reduced, and the ortho 
molecules become almost isotropically aligned. 



It is seen that the ortho and para rotational wavepack- 
ets evolve differently after the application of the pulse, 
and transient moments of simultaneous alignment and 
anti-alignment for the two species can be found. For in- 
stance, at t = 2 ps, the para molecules are aligned while 
the ortho molecules become anti-aligned. Naturally, the 
rotational energies of both ortho and para molecules are 
increased as the result of the pulse. Further selective con- 
trol of the isomers can be achieved by applying an addi- 
tional pulse (of the same intensity and duration, for sim- 
plicity) at t — 1.9 ps, when the ortho and para molecules 
are on the way of becoming anti-aligned and aligned, re- 
spectively. At this time, the alignment factors of both 
isomers are close to each other, and are a bit larger than 
1/3. The dynamics of the alignment factors after the sec- 
ond pulse is shown in Fig. 3] (again, for temperature of 
20 A^. 

It depicts a selective alignment (and anti-alignment) of 
the para molecules, while leaving the ortho modifications 
more or less isotropically distributed. Moreover, the to- 
tal rotational energy of the ortho and para molecules is 
affected in a different way by the second pulse. While 
the rotational energy of the para molecules is further in- 
creased, the energy of the ortho molecules is reduced as 
the result of the action of the second pulse. 

V. CONCLUSIONS 

We considered interaction of the water molecules with 
strong and short off-resonant laser pulses and analyzed 
dynamics of rotational wavepackets of ortho and para 
spin isomers of water. We showed that these two forms 
of the water molecules experience distinct series of tran- 
sient alignment and anti-alignment events, which leads to 
the angular separation of an initially isotropic mixture of 



nuclear spin modifications. Moreover, we showed that 
by applying a second, properly delayed laser pulse it is 
possible to suppress the rotational energy of one isomeric 
species, and leave it in a more or less isotropic state. On 
contrast, the other species will experience an enhanced 
series of transient alignments and anti-alignments. Such 
a selective alignment establishes grounds for the actual 
separation of ortho and para water isomers via interac- 
tion with additional spatially inhomogeneous static or 
time-dependent fields. This is a subject of our ongoing 
research. 

We acknowledge the support of the Israel Science Foun- 
dation. This research is made possible in part by the 
historic generosity of the Harold Perlman Family. lA is 
an incumbent of the Patricia Elman Bildner Professorial 
Chair. 



VI. APPENDIX- ROTATIONAL 
WAVEFUNCTIONS CLASSIFICATION 

The eigenstates of the water molecule are classified 
using the C2v{M) Molecular Symmetry (MS) group [H. 
Considering the water molecule eigenstate (Eq. ([1])), we 
define the nuclei rotational and vibrational, as well as 
electronic degrees of freedom as a rovibronic eigenstate. 
We denote the rovibronic eigenstate symmetry as Trve 
(within the MS group), and the nuclear spin state sym- 
metry as Tns (within the MS group). The water nuclear 
spin states are divided into either 5 = or 1 nuclear spin 
states, denoted as para or ortho, respectively. The Pauli 
principle imposes symmetry limitations on the valid 
basis functions for expressing the internal molecular 
wavefunction, ^int (the entire molecular wavefunction, 
without the translational degrees of freedom). A rovi- 
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bronic state having Trve symmetry can be only combined 
with a nuclear spin state having r„s symmetry if the 
product of the two symmetries is an allowed symmetry 
for ^int- In this Appendix we review the symmetries of 
the rotational eigenstates, and identify each rotational 
eigenstate as para or ortho. 

We start by examination of the symmetry of the 
internal wavefunction, ^mt- The character table of the 
MS group is given in Table [H The group operations are: 
E, (12), E* and (12)*, where E is the identity operator, 
(12) is the permutation of the Hydrogen nuclei operator, 
E* is the inversion operator (both for the nuclei and 
electron coordinates) and (12)* is the combination 
{12)E*. In Table U only the symmetries B2 (positive 
parity with respect to inversion) and Bi (negative 
parity with respect to inversion) change sign under the 
operation (12) (i.e. the permutation of the two Hydro- 
gen nuclei, which are fermions). Therefore, they are 
the only allowed symmetries of the internal wavefunction. 

Next, we classify the nuclear spin wavefunction. 
There are two possible spin functions for s = 1/2 proton 
spin: a = \^,^) and /3 = |^,— where a and /? 
correspond to spin up and spin down states, respectively. 
There are four possible spin configurations: 

m = 1 : aa = 

m = -1 : (3(3 = ^l^] . (17) 

One should bear in mind that the spin functions are 
invariant to E*, and also for them (12)* — (12). $1V 
and ^It] are invariant to the group operations and have 
symmetry Ai (see Table IJ). The characters of ^ns 
and (which transform together) are 2,0,2,0 (for 

E,{12),E*,{12)*, respectively), i.e. <^>if] and ^^n] gen- 
erate the representation Ai © B2. The four nuclear 
spin wavefunctions therefore generate the representation 
= 3Ai © B2. 



C2v{M) E 


(12) E* (12)* 


Rotation i?" 


r>7r r>7r r>7r 
rif, Jlc 


A, 1 


111 


A2 1 


1 -1 -1 


Bi 1 


-1 -1 1 


B2 1 


-1 1 -1 



TABLE I: The character table of the C2v{M) MS group. E 
is the identity operator, (12) is the permutation of the Hy- 
drogen nuclei operator, E* is the inversion operator (both for 
the nuclei and electron coordinates) and (12)* is the combi- 
nation (12)_E*. The effect of the MS group on the rotational 
eigenstates is equivalent to the effect of the operations of the 
molecular rotation group D2 (i.e. , , , K^) on the ro- 
tational eigenstates. 

As the next step, we classify each rotational eigen- 
state. It turns out that the effect of the MS group on 



the rotational eigenstates is equivalent to the effect of 
the operations of the molecular rotation group D2 (i.e. 
R° , Rl , , R^) on the rotational eigenstates (see Table 
m. £■ or i?'' is the identity operation. (12) (permutation 
of the two Hydrogen nuclei) is a rotation about 6 by tt 
(see Fig. [1]), which we denote as . E* is rotation about 
c by TT, i.e. Finally (12)* is a rotation about a by tt, 
i.e. The operations are given by [l|: 

R^\J,k,m) = e*'=^|J,fc,m) 

Rl\J,k,m) = (-l)-^e-2»^-"|J,-fc,m) , (18) 

where a — 0,7r/2 for b and c, respectively. Using Eq. 
(dH]) for each eigenstate, one can find the characters of 
the group operators operating on the eigenstate, com- 
pare them to Table IH and identify the symmetry of the 
eigenstate. The results are given in Table |TTT] (presented 
here only for J = 0, 1, 2 and 3). 

We assume the vibrational and electronic subsystems 
to be in the ground state (that has Ai symmetry). The 
overall rovibronic symmetry is therefore determined by 
the rotational symmetry. 

One may form a valid basis wavefunction for expressing 
^int by combining a rovibronic state having a Trve sym- 
metry and a nuclear spin state having a r„s symmetry, 
in such a way that the product of these two contains Tint 
(Trve 'S) Tns D Tint, whcrc Tint IS an allowed symmetry 
for ^int)- A detailed list of rovibronic symmetries and 
their corresponding spin symmetries is given in Table [ill 







T int 


Ai 


B2;- 


B2; Bi 


A2 


-,B2 


B2; Bi 


Bi 


-;3yli 


B2; B\ 


B2 


3A;- 


B2; Bi 



TABLE II: The symmetries of the rovibronic eigenstates ap- 
pear in the left column. The right column displays the possi- 
ble symmetries of the internal eigenstates. The middle column 
presents the corresponding nuclear spin symmetries such that 
can be either B2 or 3Ai. The '— ' sign 
indicates that no appropriate Tns exists for the combination 
of symmetries in the first and the third columns. 

From Table |TT] one can relate the Ai and A2 rovibra- 
tional wavefunctions with S = nuclei spin, i.e. para 
water (the antisymmetric 5 = wavefunction has a B2 
symmetry) . Similarly, the Bi and B2 rovibrational wave- 
functions are related with the S = 1 nuclei spin, i.e. ortho 
water (the symmetric S = 1 wavefunctions generate the 
representation 3^i). The statistical weight is therefore 
clear from the Table: there are three times more ortho 
species than para species (75% of the water is ortho and 
25% is para). 

Finally, Table IIIII identifies the eigenstates of water 
molecule as para or ortho for various values of J and 

T . 
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TABLE III: Para and ortho wavefunctions (here presented 
only till J = 3) identification. The water asymmetric rotor 
eigenstates are |J, r, m), where J is the angular momentum 
quantum number (J = 0, 1, 2, ...), and m is the quantum num- 
ber of the angular momentum projected on the laboratory 
Z axis (since there exists no preferred direction in space, we 
choose here m = without loss of generality). Index r numer- 
ates various 2J + 1 eigenstates arranged from lowest energy 
(r = —J) to highest energy (r — +J). The Table displays 
also the symmetry of the rotational eigenstates. 
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